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We study the simplicial quantum gravity in three dimensions. Motivated by the Boula- 
tov's model which generates a sum over simplicial complexes weighted with the Turaev-Viro 
invariant, we introduce boundary operators in the simplicial gravity associated to compact 
orientable surfaces. An amplitude of the boundary operator is given by a sum over trian- 
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gulations in the interior of the boundary surface. It turns out that the amplitude solves the 
Schwinger-Dyson equation even if we restrict the topology in the interior of the surface, as 
far as the surface is non-degenerate. We propose a set of factorization conditions on the 
amplitudes which singles out a solution associated to triangulations of S 3 . 
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1. Introduction 

In order to understand a quantum theory of gravity, we need to clarify what we mean 
by integrations over Riemannian metrics on a <i-dimensional differentiable manifold M. 
One of the possibilities proposed by Regge Jl| is to discretize the integrals and and replace 
them by sums over smooth triangulations on M. In this approach, each triangulation T 
specifies a Riemannian metric on M in the following way. We impose that an interior 
of each <i-simplex in T is flat Euclidean and that the curvature has a support only on 
(d — 2)-simplexes. Moreover all the edges in T are regarded as being straight and having 
the same geodesies length a. (There is another version of the theory where the edges are 
allowed to have different lengths 0, but we do not examine this possibility here.) These 
conditions are enough to define a Riemannian metric associated to the triangulation T of 
M. It is reasonable to expect that any metric on M can be approximated in this way 
by using sufficiently large number of simplexes and by taking the size a of each simplex 
to be infinitesimally small. It is then important to know what types of triangulations 
give relevant contributions to the summations. If smooth triangulations with infinitely 
many simplexes dominate the summations, it would make sense to define the integrals 
over metrics by the sums over triangulations. In two dimensions, this program has worked 
rather well and we have learned quite a lot about two-dimensional quantum gravities from 
this approach for the last few years. The purpose of this paper is to study the three- 
dimensional case. 

Let us be more specific about what we mean by the sums over triangulations. In 
order to quantize the Einstein gravity, we would like to perform integrals over metrics g^ 
weighted with the exponential of the action given by S[g] = A J d d Xy/g + -^ J d d x^fg~R 
where R is the scalar curvature. Suppose that the metric g is associated to some trian- 
gulation T of M. Since each <i-simplex has the same volume proportional to a d , the total 
volume J d d Xyfg of M should be proportional to a d rid(T) where rii(T) (i = 0, 1, ...,d) is 
the number of i-simplexes in T. The integral of the scalar curvature R on the other hand is 
expressed in the cases of d = 2 and 3 as follows. When d = 2, the curvature has a support 
on the vertices of T. Since the deficit angle around each vertex v is given by (2-7T — ^t^(v)) 
where ta(v) is the number of triangles containing v, we obtain 

fd 2 x^g-R~ J^ (l-±r A (v))=no(T)-±n 2 (T). 

v.vertices ^ ' 



Especially when the manifold M is closed, 3ri2(T) = 2n\(T) since every edge on T is 
shared by exactly two triangles. In this case, the right-hand side in the above becomes 
X = no(T) — ri\{T) + ri2(T), the Euler number of M. When d = 3, the curvature is 
concentrated on the edges of T and is characterized by the deficit angle around each edge 
e given by (2tt — cos _1 (^)rt(e)) where Tt(e) denotes the number of tetrahedra containing 
e. Since each edge has the geodesic length a, the integral of the scalar curvature is given 
by 



/ d^x^/gR ~ a } 
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2nn 1 (T)-6cos- 1 [ - ) n 3 (T) 



Therefore, for both d = 2 and 3, the action S[g] is expressed as a linear combination of 
n d (T) and n d _ 2 (T). 

S[g] = -Kn d -2{T) + (3n d (T). 

The partition function of the simplicial quantum gravity is then defined by 

1 



E 



C(T) 



D Kn d - 2 (T)-(3n d (T) 



where the sum is over all possible triangulations of M, and the factor C(T) is the order of 
a symmetry of T if any. This factor is included in order to take into account the volume 
factor of the diffeomorphism group. First of all, we need to know if this summation is 
convergent at all. Since Z d can be written as 



E Z n d ,n d - 2 (M)e 

rid,n d _2=0 



K,n d ^ 2 -0ri d 



where Z. 



n d ,n d -2 



(M) is a number of triangulations of M with given values of n d and 
n d -2 weighted with the factor C(T), we would like to know the asymptotic behavior 
of Z nd7 n d _ 2 (M) for large values of n d and n d -2- In two dimensions, n<i and no are not in- 
dependent and they are related as as no = \ri2 +X- Especially when x = 2 (i.e. M ~ S 2 ), 
the exact asymptotic behavior of Z n2 ^ no is known from the combinatorial study of ref. 
as 

Zn„.n (5 a )~wJ'e /3cna (n 2 
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J n 2 ,n 
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oo, n = -n 2 + 2) 
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The sum over triangulations can also be evaluated by the matrix model [ffl, and the above 
formula for S 2 is reproduced by the sum over one-particle irreducible vacuum diagrams 
H|] in the large- N limit of the matrix model || || @. (The sum over connected vacuum 
diagrams also gives a similar formula with a different value of (3 C ; e^ c = 2 • 3 3 .) When M is 
a surface of higher genus (x < 0), subleading corrections || to the large- N matrix model 
give 

Z n2no (M) ~ n~ ix ~ l e^ n * (n 2 -> oo, n = -n 2 + X ). 

Therefore the contribution to Zd=2 from T with many triangles is expressed as 

Z d = 2 ~e*xJ2 n 2 iX ~ le ~ ~ Pc)n2 ~ 

where (3 = (3 — |«, and it is convergent for J3 > j3 c . Moreover as (3 approaches to the 
critical value (3 C , the contributions from triangulations with n 2 >> 1 become more and 
more significant for Z ( j i=2 with x < and for d%Zd=2 with % = 2. Thus we have a good 
reason to believe that the program of the simplicial quantum gravity works well in two 
dimensions. Indeed this approach has turned out to be quite successful and its validity 
has been confirmed by comparisons of the matrix model computations with the results of 
the continuum field theories [0 |TT[ @ |TJ] |TJ QT5| |T§. 




Recently several numerical results have been reported on the three-dimensional sim- 
plicial quantum gravity |17j |H| [|19| [^0[ plU - Their results seem to be compatible with 



the following bound on the number of triangulations of S 3 . 

Zn 3 (S 3 ) =^4 3ni (S 3 ) < Ce^ (n 3 - oo), (3) 

Til 

for some constants C and (3 C . In three dimensions, the Euler number \ — n§—n i+n2— n3 is 
zero for any triangulation of a closed manifold. We also have n 2 = 2ns which follows from 
the fact that any triangle in T is shared by exactly two tetrahedra. Therefore 77,3 = ni — no, 
and 77.3 should always be less than n\. On the other hand, since each edge in T belongs 
to at least one tetrahedra, n\ must be less than or equal to 677.3. Namely 773 < n\ < 6773. 
Thus, if the exponential bound (|3|) is indeed correct, the sum over triangulations in eq. (TJ) 
would be bounded as 



^=3 < Yl Z nzni e-^-^ n * < Cj2e~ 



^ -(/3-6K-/3 c )n 3 

"3, Til «3 



for k > and 

Z d=3 < J2 Z n3ni e-^-^ < Cj2e~^-^ 

n 3 ,n 1 n 3 

for k < 0. The sum ([!]) would then be convergent at least for j3 — 6k > (3 C (k > 0) and 
for (3 — K > P c (0 > k). It is thus important to develop an analytical tool to examine the 
behavior of Z na (S 3 ). 

In two dimensions, the matrix model [ffl has been found useful to study the simplicial 
quantum gravity, and there have been some attempts to extend it to three dimensions 



||22|| [23] 1 24]. In this so-called tensor model, one considers and integral over rank-three 
tensors rather than matrices, and the perturbative expansion of the integral generates 
simplicial complexes, i.e. collections of tetrahedra whose boundary triangles are pairwisely 
identified. However the collections of tetrahedra generated in this model are not necessarily 
triangulations of smooth manifolds. In fact their Euler numbers are not always zero, but 



in general non-negative integers [22 



The basic reason for the failure of the tensor model is that the model does not contain 
as many parameters as one needs to distinguish various topologies in three dimensions. In 
two dimensions, a topology of a closed orientable surface is characterized by the number 
of handles, and two parameters in the model are enough to control both the topology 
and the volume ni of the surface. The matrix model indeed has two parameters, the 
size iV of the matrix and the coupling constant g. In three dimensions, on the other 
hand, we need at least three parameters to control three independent variables 713, ni and 
X = no — n\ + ri2 — n^. The simplest version of the tensor model has only two parameters, 
the size iV of the tensor and the coupling constant, and it is not possible to force x to be 
zero. 



Recently, Boulatov proposed an improved version of the tensor model |25[] by incor- 
porating a gauge symmetry of a finite group or a quantum group. If one takes the gauge 
group to be a finite group G, the order of G couples to the Euler number of the simplicial 
complex. For example, when G is the cyclic group Z p , the partition function of the model 
is expressed as a sum over simplicial complexes T weighted with a factor p x ( T '\Hi(T, Z p )\, 
where \Hi(T; Z p )\ is the order of the first cohomology group of T with coefficients in Z p . 
In three dimensions, the necessary and the sufficient condition for a simplicial complex T 
to be a manifold is x(T) = 0, and otherwise x(T) > 0. Thus, if one could make sense of the 
limit of p — > 0, the perturbative expansion of the model in this limit would be dominated 
by the sum over simplicial manifolds with x(T) = 0. 



The p — > limit of the Z p model not only imposes x = but gives an additional 
condition on topologies of simplicial manifolds, i.e. the first and the second homologies 
of the manifolds are restricted to be trivial [E5| . Without any condition on the topology, 



the number of three-dimensional simplicial manifolds with a given number 77,3 of tetrahe- 
dra grows at least factorially in 77,3 |22| , and the sum over simplicial manifolds becomes 
divergent. Thus one needs to impose restrictions on the topologies. It is not clear if the 
conditions on the homologies are enough to tame the factorial growth since there are many 
topologies in three dimensions with the same homologies as those of S 3 . If not, the sum 
over triangulations generated by the Z p model would still be divergent, and we would need 
to develop a mechanism to impose further restrictions on the topologies. 

In this paper, we examine the Boulatov's model associated to the quantum group 
SU q (2). In this model, the sum over simplicial complexes is weighted with Aq I q (T) 
where 

2(fc + 2) 2ni/(k+2) 

^1 ( g l/2 _ g -l/2)2 ' y 



and I q (T) is the Turaev-Viro invariant [ 26fl . We define boundary operators of the model 
associated to compact orientable triangulated two-dimensional surfaces, and derive the 
Schwinger-Dyson equation for the amplitudes of the operators. It turns out that the 
Schwinger-Dyson equation takes the same form for any value of q when the surfaces involved 
are non-degenerate. Since the amplitude of the boundary operator in the Boulatov model 
is expressed a sum over triangulations in the interior of the boundary surface weighted with 
A^I q depending on q, the independence of the Schwinger-Dyson equation on q suggests 
that we can restrict the topology of T in the interior of the surface without spoiling 
the equation. Indeed we find that, associated to an arbitrary closed orientable three- 
dimensional manifold, one can construct a solution to the equation. In order to single out 
a solution associated to triangulations of S 3 , we introduce a set of factorization conditions 
for degenerate surfaces with which the equation can be solved uniquely and inductively. 

In the case of the two-dimensional simplicial gravity, the Schwinger-Dyson equation 
can be derived from purely combinatorial considerations without using the matrix model. 
This is also true in the three-dimensional model studied here. However, in order to take 
into account the symmetry factor to each graph, it is more convenient and transparent to 
employ the perturbative expansion of the Boulatov model. 

This paper is organized as follows. In Section 2, we briefly describe of the Schwinger- 
Dyson equation in the large-iV matrix model and its combinatorial meaning. We then 



define the Boulatov's model in Section 3 and examine its properties. The boundary oper- 
ator are defined in Section 4, and the Schwinger-Dyson equation is derived. In Section 5, 
we show that there is a solution to the equation associated to each closed orientable three- 
dimensional manifold. The factorization conditions on the amplitude which characterize 
the solution associated to S 3 are defined in Section 6. The Schwinger-Dyson equation 
combined with the factorization conditions is shown to have a unique solution. 



2. Schwinger-Dyson Equation in the Matrix Model 

Let (Mij) be a N x N hermitian matrix and consider the following integral. 

Z= f[dM]exp{-S{M)) (4) 

where 

S(M) = -trM 2 ^=trM 3 . 

by expanding the integral ([|) in powers of g, Z is expressed as a sum over orientable closed 
trivalent graphs, each of which is dual to a triangulation of an orientable closed (but not 
necessarily connected) surface. A sum over closed and connected simplicial manifolds is 
obtained by taking the logarithm of Z as 

logZ = V — — ( 9 ) jW T ) = V -!—n^( T )N^ T \ 

Especially the sum over triangulations of S 2 is given by the leading term in the large- N 
limit of log Z. 

The basic quantity we examine in this section is the expectation value of trM n at the 
large- N. 



Un N 



Jim -J-- (trM n ) = lim -^ ' 4 f[dM]trM n exp(-S(M))- 



The power of N is multiplied to (trM n ) so that u n gives a finite value in the limit of 
iV — > oo. The Schwinger-Dyson equation for u n [^7[, [^SJ is derived from 



= ^j\dM]^- [(M- 1 ), j exp(-5(M))] 



n-2 (5) 

= -(trM n ) + -j=(trM n+1 ) + ^(trM k trM n - k - 2 ). 

v^ k=Q 
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Here the integration by part is justified in each term in the perturbative expansion in g. 
In the limit of large- N, the last term in the right-hand side of eq. (|5p factorizes as 

since only the planar graphs contribute. Thus we obtain a recurrence relation for u n as 

n-2 



gU n+ i + 22 U kU n -k-2- (6) 



u n 

fc=0 

This is the Schwinger-Dyson equation in the large- N limit of the matrix model. 

This equation has the following simple graphical interpretation. The expectation value 
u n can be expressed as a sum over planar trivalent graphs with n external legs. The graphs 
are not necessarily connected, but they do not contain a closed disconnected part because 
of the factor 1/Z in u n . Let us take one of the graphs T in the summation and pay attention 
to one of the external legs on it. The leg should be either connected to a trivalent vertex 
in T or to another external leg. In the former case, the rest will be a graph with (n + 1) 
external legs. If on the other hand the leg is connected to another external legs, the line 
connecting these two legs will separate the original planar graph into two disconnected 
graphs with k and (n — k — 2) external legs (k = 1, 2, ...,n — 2). The two terms in the 
right-hand side of the Schwinger-Dyson equation correspond to these two possibilities. 

It is also straightforward to show that the Schwinger-Dyson equation completely de- 
termines the power series expansion of u n in g. By substituting 

oo 

u n(g) = ^2 u »,»afl ,na 

712=0 

into eq. (|6|), we obtain a recurrence relation with respect to n and n^ as 

n — 2 ri2 

Un,n 2 — u n+l,n 2 -l + /_^ 2_^, U k,sUn-k-2,n 2 -s- (7) 

fc=0 s=0 

Since uo(g) is normalized to be 1, the initial condition of the recurrence relation is uq j712 = 
£0,7125 with which u n ^ n2 is determined inductively by solving (|7|) following the bibliographic 
ordering: (n, n<i) >- {n\ n' 2 ) if n 2 > n' 2 or jt-2 = n' 2 , n > n'. 

Since the nonlinear term in the Schwinger-Dyson equation (|6|) is the convolution of 
ii n 's, it can be transformed into the integral equation 



where 



u r 



dXp(X)X r 



(9) 



and P J denotes the principal part of the integral. This equation can be solved using the 
technique of the Riemann-Hilbert correspondence || as 

p(\) = l[l + lg(a + b)+g\W(\-a)(b-\) (a < A < b) 

TV I 







where a and b are determined by 



(A < a or b < A) 



g(b-a) 2 + 2(a + b)[2 + g(a + b)} =0 , (b - a) 2 [l + g(a + b)} = 1. 

We can also compute the partition function (|]) itself from the solution to the 
Schwinger-Dyson equation. By using the relation between u% and Z, 



u 3 



1 . , -. 3 1 dZ 

—(trM 3 ) = ^r — — 



(10) 



we obtain 



n 2 r 9 

Z(^) = Z(0)exp(— / dg'us(g')). 
<j Jo 



The approach described here using the Schwinger-Dyson equation is in fact equivalent at 
the large-iV to the standard computation || using the eigenvalues Xi (i = 1, ...,N) of 
M. By inverting the integral transformation (^|), we can express p(A) in terms of u n as 

P(A) = /| e - ± -1„„ A » = / | e - ( ^ rexp(8fc il )> 



T1 = 



- 1 



A 7 



))• 



Namely p(A) is the density of the eigenvalues of M. 



3. Boulatov's Model 



In this section, we describe the model due to Boulatov |25|] and examine its properties. 
Let us consider the following integral 



[dM] exp(-5(M)) 



11) 



1 

s ( M ) = a E E E i M 4^r 2a3 i 2 

H,a 2 ,a3=l {ii,j 2 ,j3} {-Jfc<w.fc<ifc} 



with 5(M) given by 

AT 

6 

a 

_^i V V V r-DEL.iJ^ -? 2 M (12) 

ai,...,a 6 =l { 3l ,...,j 6 } {-j k <m k <j k } v '9 

v (— tf)5-/fe=l mfe ]\/f3lhhi a l a 2a,3 T\ / rJ3J4,J5;a-3a4.a-5 JJfh 35 J6;«l«5«6 j^ / rJ2J6J4,;a 2 a 6 a 4 , 
\ "i — mi— 7712— 7773 7773 — 77747775 7771—77757713 7712 — 77767774 

where 

' Jl J2 J3 

J4 h k 
is the Racah-Wigner 6j-symbol of the quantum group SU q (2) [[29] with q = e 2 W( fc + 2 ) anc j 

= 2(fc + 2) 

9 ( g l/2_ g -l/2)2- 

The integration variable M^^ 3 ^ a2a3 obeys the conditions 

(1) M™)™Tata 2 a, = unless \ 3l - j 2 \ < j 3 < j x + j 2 and 3l + J2 +j 3 <k 

(9\ -Tj3\hJ3\a 1 a 2 a z _ . \j x +j 2 +j 3 , _ x r n 1 +m 2 +m 3 M Ji j 2 jz\^a 2 a z 
\ A ) lvl m 1 m 2 m 3 ~ \ L ) \ H) JW -mi-m 2 m 3 

CX\ M Jlj2j3 ' aia2a3 = M j2j3jl,a2a3ai 

\"/ 777177727773 777277737771 

The perturbative expansion of this integral gives a sum over Feynman graphs each of 
which is dual to an oriented simplicial complex in three dimensions. A contribution of 
each simplicial complex T is evaluated by assigning a triplet of indices (a e ,j e , m e ) to each 
edge e in T and by summing over the indices on all the edges with the conditions 

lid (A) -ie 2 (A)| < Je 3 (A) < Jei(A) + Je 2 (A) , J ei (A) + Je 2 (A) + Je 3 (A) < k 

for every triangle A in T where e r (A) (r = 1, 2, 3) are the edges of A. We then obtain 

^z q =Y,^9Kr {T) ^ {T) K o{T) - 1 ^n 

where the sum J2t * s over connected simplicial complexes, and I q (T) is given by 



w=T^n£ n (-i) 2 ^ e +i] 9 

P {j} e:edges 

,.^LL^ ' \je*(t) je 5 (t) je a (t) J 



where 



gi+1/2 _ q-J- 1 / 2 

i 2 J + ^9 = g l/2_ g -l/2 



and e r (£) (r = 1, ..., 6) are the six edges of the tetrahedron £. 

It has been shown by Turaev and Viro |26j that I q (T) is a topological invariant of 
T, i.e. IqiTx) = IqiTz) if ^1 an d T2 are combinatorially equivalent. Especially when T 
is a triangulation of an orientable manifold M, it has been shown in |3(J |31]] that I q (T) 
is equal to |r g (M)| 2 where r q {M) is the Witten-Reshetikhin- Turaev invariant |5!|. (The 



relation between I q {T) and r g (M) at q — » 1 has also been studied in [|33| ]34[]. See also [J35] 
and J36|.) For example 

7 g (S 3 ) = l, I^S 1 x S 2 ) = A q . 



In general, it has been shown by Kohno [B7| that I q (M) is bounded from above as 



7 g (M) < Aj (M) (14) 

where h(M) is the Heegaard genus. Any three-dimensional manifold can be decomposed 
into a pair of handlebodies Tin of genus h glued together on their boundaries (the Heegaard 
splitting), and the Heegaard genus h{M) is the minimum number of such h. Especially 
h(M) = if and only if M is homeomorphic to S 3 and otherwise h{M) > 1. 

By introducing a new variable fiby N = /i/A g , the partition function Z q is expressed 
as 

lo §^ ^E^^^WT). (15) 

Suppose that we could take the limit q — ► so that A g — ► while maintaining the inequality 
(P^[). As we mentioned in the introduction, x(T) > for any simplicial complex T in three 
dimensions and the equality holds if and only if T is a simplicial manifold. Because of the 
inequality ([T3]) , the q — > limit of log Z q would be dominated by a sum over simplicial 
manifolds with x(T) = and h(T) = 0, i.e. T ~ S 3 . Unfortunately, it is not clear how 
to make sense of this limit since the model is defined only when q is a root of unity. In 
the following, we will develop an alternative mechanism to impose the condition T ~ S 3 
without using the q — > limit. 



4. Schwinger-Dyson Equation in Three Dimensions 

10 



In the case of the matrix model, the integral over the hermitian matrix M in (|l]) can 
be reduced to an integral over the eigenvalues Xi (i — 1, ..., N) of M 

N 

J[dM] eMS(M)) = J J] d\ t \\{K - A,) 2 exp(- J^A 2 - ^=Af )), 



3VN 

and the powerful methods of the WKB approximation [f3J and the orthogonal polynomials 
have been used to successfully analyze the model. This approach, however, does not 
seem to have an obvious generalization to three dimensions. It is not clear if the integral 
over M^^^ a2 ° 3 could be significantly reduced by introducing some variables similar to 
the eigenvalues. 

In Section 2, we saw that the large-iV matrix model can also be studied by using 
the Schwinger-Dyson equation for u n , the expectation value of trM n . In fact u n is the 
integral transform @ of p(\), the density of the eigenvalue. This is not too surprizing 
since both the WKB approximation and the Schwinger-Dyson equation at large- N make 
use of the classical equation of motion of the model. The point here is that the method of 
the Schwinger-Dyson equation is applicable to the three-dimensional case. 

In case of the matrix model, the expectation value of trM n is expressed as a sum 
over planar graphs with n external legs. Thus trM n may be regarded as an operator 
which creates a circle with n edges in each of the planar graphs. The natural thing to 
consider in three dimensions should then be an operator which creates a closed connected 
triangulated surface E. Indeed we can construct such an operator as follows. We attach a 
triplet of indices (j e , m e , a e ) to each edge e on the surface S and associate M-jJ-^|^ aaa3 
(jr — ie r (A) 5 "V = ± ?ri e r (A) 5 a r = &e r (A)) to each triangle A. The signs (±) of m r for a 
pair of triangle sharing the edge e are chosen to be opposite. We then define a polynomial 
F-z(j; M) of M by summing over {m} and {a} as 

F E (j;M)= J2 II ^ 1 [(-?r e + (-?)" me ] 

{m},{a} eiedges /-.^s. 

X IT yiy-ieiCA)Je 2 (A)ie3(A)a ei (A)a e2 (A)ae3(A) 

11 ±m ei (A)±™ e2 (A)± r ". e3 ( A) 

A:triangles 

By taking an expectation value of F%, we obtain a function ^t,(j) of the indices {j} on E 



as 



*s(i) = (Fs(i; M)) = i f[dM]Fv(j; M) exp(-5(M)). (17) 



11 



We have not yet summed over the indices {j}. 

We saw in eq.(|13D that the weight I q (T) for the simplicial complex T in in the summa- 
tion (|15j) is given by a sum over {j} on the edges of T. We may regard I q (T) as a partition 
function of a lattice statistical model on T. In such a topological lattice model, we can 
introduce a notion of physical states on the closed triangulated surface E as follows |26j [j33 | 



34}] . We start with a space Ve spanned by functions $(j) of {j} on E and call it a space of 
states. In order to define physical states in Ve, we consider a three-dimensional simplicial 
complex Tj^o s x which is homeomorphic to E x [0, 1] and whose boundaries E x {0} and 
E x {1} are triangulated as E and Ei. We then define a map -Pe ,Ei from Ve to V^ 1 by 

fto,Ex = <&(i) e Vfc - ^$(/)Pe ,e 1 (/,j) e V^ 

O'} 

where Pso.EiO'ji) is defined by attaching an index J e on each edge e on Te 0) Ei and by 
summing over {J} as 

i^o,Si(j / ,j)= noko,^) E II *( J «»i«) 

A 9 {J} e6Sx{0} 



X 

eeEx 



n ( (-d^P+iJ n (-i)"-w +u. us) 

Ex{l} VV ; L Je lq/ e£T So , El 

rr (_i)ELi J = fe (o/ J ^iw J ^w J - 3 (t) \ 



X 

t€T S() , El 



Because of the topological invariance, -Pe ,Ei is independent of the choice of the triangu- 



lation in the interior of Te ,Ei and satisfies |26 



J2 P XoMJ,3') P ^,^ti',3") = Peo&UJ") (19) 

{j'} 

where Eq, Ei and E2 are all homeomorphic to E. Especially when Eq = Ei = E, Pe,e acts 
as a projection operator on Ve- The physical subspace V^f ys is defined as a subspace of 
Ve projected out by Pe,e as 

if} 

Because of the property (|T9|), the operator Pe ,Ei defines an invertible map from a physical 
state in Ve 2 onto a physical state in Ve - Therefore V^ ys and V^f are isomorphic. 
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Since ^s(i) given by ( |TTD is also a function of {j}, we may consider its pairing with 
a physical state $ on E as 

($,^)=A| n2(E)+,l(E) $:$(j)^0-) (20) 

where /t(E) is the genus of E. The factor A| 2 will simplify the equations in the 

following. As we shall see below, $e(j) itself is not a physical state, and the pairing 
($, ^s) depends on the triangulation E of the surface. The Schwinger-Dyson equation 
describes the dependence of ($, ^s) on E. 

Among physical states on E, especially important are the Hartle-Hawking type states 



34|1 defined as follows. A handlebody Tin of genus h is obtained by cutting out 2h pairs of 
discs from a boundary of a three-dimensional ball B 3 and by identifying them pairwisely in 
opposite orientations. The boundary of Tih is then a closed orientable surface of genus h. 
The boundaries of the discs become simple homology cycles on the boundary of Tih which 
do not intersect with each other and which are contractible in Tih- They are called the 
meridians of the handlebody. Let us consider a simplicial manifold n^ homeomorphic to 
a handlebody of genus h whose boundary is triangulated as on E (/i(E) = h) and whose 
meridians are given by {ax, ...,cth} where each on is chosen to be a sequence of edges on 
E. The Hartle-Hawking type state <&( a >(j) associated to the handlebody is defined by 

K s V}^ Vl ' L J iq/ eew<r> (21) 



hh. ^ I J e 4 (t) J e 5 (t) Je 6 (t) J ( 



X 

ten({ a }) 

Because of the topological invariance, $w(j) defined in this way is independent of the 
choice of the triangulation of the interior of the handlebody, and it gives a physical state 
on E. 

By expanding integral ( (FTP in powers of g, the pairing of \&e with the Hartle-Hawking 
type state $( a ) is expressed as 

{j} 



ni(T) 



xA^^/rfu^-T)) 

= £ (§P)^ (T) A* Bl(T) Aj (T) -* x(E) i ff (wWu E (-r)) 



(22) 



T:9T=E 
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where the sum J^ T is over simplicial complexes whose boundary is E and I q (Hj? U^ (— T)) 
is the Turaev-Viro invariant (|13|) for the simplicial complex obtained by gluing n^ and 
T together on their boundaries after inverting the orientation of T . As in the case of the 
matrix model, the dual graph of T is not necessarily connected though it does not contain 
a closed disconnected part. Especially a triangle on the boundary E of T may be attached 
to another triangle on E. For any simplicial complex T bounded by E, its Euler number 
x(T) is greater than or equal to |%(E), and the equality holds if and only if Hj? U (— T) 
is a simplicial manifold. If we could take the limit q — > 0, the right-hand side in the above 
would become a sum over T such that x(T) = |x(E) and W^ Ue (— T) ~ S 3 . Although 
this limit itself is not well-defined, we will find a set of factorization conditions which 
together with the Schwinger-Dyson equation imposes the similar conditions on T. 

In order to derive the Schwinger-Dyson equation, we choose a triangle Aq on E and 
consider an open surface Ea obtained by removing the triangle Aq from E. Associated 
to E Ao , we define a boundary operator F^ Aq (j; M)%**?* ma by 

3 

F S a;M)^X 3 m 3 = J2 H S ( m r,m eAAo) )S(a r ,a eAAo) ) 

{ m },{a}r=l 

x II { - z y 1 [(-s) me + H?r me ] 

eiedges 
X If 7lrf-Jei(A)Je 3 (A)Je3(A);a el (A)O e3 (A)Oe3(A) 

11 ±m ei ( A )±m e2 ( A )±m e3 (A) 

^■.triangles 

As in the case of the matrix model, we make use of the identity 

\\^ E dM nLa^a s IF^M)T^*M-S{M))] =0 

^ r l r 1 'JlVlm^m'ym-x 



{m},{a} 



where we have set j er (A ) = 3r ( r = 1) 2, 3). By performing the derivative with respect to 
M, we obtain 

. Z — '. I J4 35 36) ^o 

34,35,36 

<9Fv (j- M) aia2a3 ^ ' 

. V^ / __^0 rLl fT r n 1 ra 2 rao, , 
/ -> ^ Q/\ / /-jiJ2j3;«ia 2 a3 ' 

{m},{a} 01Vl mi m 2 m 3 

where E^ is a triangulated surface obtained from E by decomposing Ao into three small 
triangles, and J4,,jb,je are spins on the three additional edges on E^ . If it were not for 
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the last term in the right-hand side of eq.(^3|), vI/e would transform like a physical state 
under the change of triangulation E — > E^ . 

Since each M in F-£ A (j;M) corresponds to a triangle A on Ea , the last term 
(-^■F-£ A (j; M)) in the right-hand side of (^) can be expressed a sum over closed surfaces 
Sa,a obtained by removing A and Ao from E and by identifying their edges. Let us 
denote the number of common edges of the two triangles by E(A, Ao) and the number of 
common vertices by V(A, Ao). We can classify the position of A on Ea according to the 
values of E and V as follows. 

(1) When E(A, A ) # and V(A, A ) = E(A, A ) + 1, A is in the neighborhood of A 
on E. In this case, Ea,a has the same topology as E. 

(2) When E(A,Aq) = and V(A, Ao) = 0, 1, the removal of the triangles will create a 
new handle on E and the edges of Ao becomes a homology cycle around the handle. 

(3) When V(A, Ao) > E(A, Ao) + 2, the triangles A and Ao wrap around a cycle on E, 
and Sa,a is degenerate around the cycle. 

Let us consider a generic situation when E(A, Ao) + 1 > V(A, Ao) for any triangle A on 
Ea , which includes the cases (1) and (2). We also assume here that the surface E is not 
degenerate in the neighborhood of Ao- The cases involving degenerate surfaces will be 
examined in the next section. 

In the generic situation, by multiplying an arbitrary physical state $ on both sides of 
eq.(^3D and by summing over {j}, we obtain 

(d>, *e) = g(*, *e 2 ) + Yl /* JS(A ' Ao) (*»*Ea,a ) 

A " A;E(A,A )/0 

+ Yl (^[ $ ]a,A ,^E a , Ao ). 

A;E(A,A )=0 

The second term in the right-hand side is a sum over A in the neighborhood of Ao, and 
Sa,a is °f the same topology as E. On the other hand, /i(Ea,a ) — h{E) + 1 in the last 
term, and «a,a is defined by 

i[$}A,AoU) = J2 ( II 5 0'e r (A), ie r (A )) 1 ®(j)- 

{j e ,.(A)}f, =1 Vr=l / 

It is easy to show that this gives a map from V-f into V^ 1 ys . 
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Especially when $ is the Hartle-Hawking type state $( a ), i[$^]a,a i s a ^ so the 
Hartle-Hawking type state associated to a handlebody whose boundary is Sa,a an d whose 
meridians are {«i, ..., a^), dA } 

^[$ (Q) ]AAoO') = * (QU{9Ao}) 0') (25) 

where <9Aq is the boundary of Aq. The equation (|24]) for $ = $( a ) can then be written as 

<E<«>)=s<Eg>+ E ^ (A ' Ao) (sk Q) Ao )+ E ( S £A { o aA ° }) ) ( 26 ) 

A;E(AAo)#0 A;E(AAo)=0 

where 

(E («) } = (*(«), * E ) = A| n2(E)+/t(E) ^$(«)(j)(F E (j;M)). (27) 

{i} 

We call this the Schwinger-Dyson equation for the amplitude (E^) of the triangulated 
surface £ with the choice {cti} of simple homology cycles. We see that, when £a,a is 
non-degenerate, the Schwinger-Dyson equation (|2^) does not depend on the value of q 
explicitly. 

5. Restriction on Topologies 
5-1. General Case 



We have found that the Schwinger-Dyson equation (26) takes the same form for any 
value of q as far as the surface £ in the equation is non-degenerate. On the other hand, 
the amplitude (S^') which solves the equation is expanded as a sum over T weighted 
with A q 2 IqiTtj?' Us (— T)) as in eq. (|2~2|). Depending on the topology of T, the 
weight takes various functional forms in q. Since the Schwinger-Dyson equation for non- 
degenerate surfaces is linear in (£(")), a sum over a restricted class of triangulations T 
such that the weight takes the same functional form should also give a solution to the 
equation. 

Here we claim a stronger statement that, for each choice of a compact orientable 
three-dimensional manifold M, a sum over triangulations T such that 

?4 a) U E (-T) ~ M (28) 
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satisfies the Schwinger-Dyson equation. Being a triangulation of the manifold, T should 
also satisfy %(T) = |%(E). Let us define 



! sW >i« = £ (g§) 9 "' (T V«"' (T » 



T:?4 a) U E (-T)~M 



= X] Z ^^i(^ a )\M9 na ^ 7 



(29) 



,"3 ,,"1 
n 3 ,ni 



where the sum Ex is over triangulations T satisfying the condition (P5|), and Z„ 3ini (E; a)| M 
is the number of such triangulations with n^(T) = n^, ri\(T) —n\. We are going to show 
that this gives a solution to the Schwinger-Dyson equation. 

By definition, for each triangulation T counted in the above, n^ Ue (— T) give a 
triangulation of M. Let us pay attention to one of the triangle Aq on the boundary £ of 
T. In the simplicial complex T, the triangle Ao may be either attached to a tetrahedron 
or identified with another triangle A on £a - 

If Ao is attached to a tetrahedron in T, we can remove the tetrahedron from T and 
give it to the handlebody n^\ without changing the triangulation of M. After removing 
the tetrahedron, T becomes a triangulation T' consisting of (ns — 1) tetrahedron whose 
boundary is E-r . Obviously T' obeys x(T') = |x(E^ ) if%(T) = i%(E). The handlebody 

n£ becomes 'H^i after adding the tetrahedron on it. Since the triangulation of M 

2 

remains unchanged, T' should satisfy 

ft£2 U E ^ (-T') ~ M. (30) 



A, 







Such T" is counted in Z n3 _i ini (S^ ; ck)|m- Conversely, from any triangulation T" satisfying 
the above, we can construct a triangulation T such that the condition (|28|) holds. However 
the correspondence between T and T' is not necessarily one to one when T and T" have non- 
trivial symmetries. We shall see below that this is taken care of by the factors C(E)/C(T) 
and C(Z~ o )/C(T'). 

Let G(S) be the symmetry of E under interchange of the triangles on it, and G(Ea ) 
be the invariant subgroup of G(E) which keeps the position of the triangle Aq fixed. 
Then [G(E) : G(Ea )] = C(E)/C(Ea ) gives the number of the triangles on E which are 
conjugate to A under the action of G(E). Suppose that mi(T) of them are attached to 
tetrahedra in the interior of T, and m,2(T) of them are attached to triangles on Ea so 
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that rai(T) + m 2 (T) = C(E)/C(E Ao ). The contribution C(E)/C(T) of the triangulation 
T to Z n3ni (E; a)|M can then be divided into two parts as 






Let us show that the sum 



CPW-E(W) (3D 



T 



VC(T) 



over triangulations T with satisfy eq. (|28p and n^T) = 77,3, ni(T) = ni is equal to 
■^n 3 -i,ni(E^ ;«)|m- If we remove a tetrahedron attached to one of the mi(T) triangles 
conjugate to Ao, we obtain a triangulation T" consisting of (77.3 — 1) tetrahedra. It may 
happen that there is more than one ways to remove a tetrahedron from T to obtain the 
same triangulation T'. Let N(T — > T') be the number of ways to obtain T" from T. 
Obviously m\(T) — Y^t' N(T — ► T"), and the summation given by eq. (|31| ) is equal to 



^)-EE(^p). 



T T> 

The number N(T —> T') can be expressed as 

C(T) 



AT(T -»■ T') 



C(T, T') 

where C(T, T') is the order of the symmetry of T which keeps the position of the tetrahe- 
dron at T \ T' fixed. Therefore eq. ( |5T] ) can be expressed as 

On the other hand, 

/C(££ o )/C(£ Ao ; 
Zng_i jni (E^ o ;a)| M = C(E Ao ) ■ 2^ I £77^7; 

where N(T *— T") is the number of ways to construct T from T" by attaching a tetrahedron 
on the boundary of T' . The number N(T <— T") can also be expressed as 

C(T') 



N(T' <- T) 



C{T, T') 
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ining these, we obtain 




^n 3 -i,ni(S^ o ;a)| M 


= C(E Ao 




= ^(Sa 



\f\c{T,T>) 

V f m 1 (T) \ (32) 

V v c(r) ; ' 

Next we consider the case when Ao on the boundary of T is identified with another 
triangle A on Ea - If A is in the neighborhood of Ao on E, i.e. i?(A,Ao) 7^ and 
1/(A, Ao) = i?(A,Ao) + 1, we can remove the triangle A = Ao from T to obtain T" 
whose boundary is £a,a • I* 1 this case, the topologies of the boundaries are the same, 
i.e. x(^) = x(Sa,a )- Since T' is obtained from T by removing one triangle, E(A, A ) 
edges and (V(A, Aq) — 2) vertices, the Euler number of T' remains the same as that of T; 
X(T') = X (T) - 1 + E(A,A ) - (V(A, A ) - 2) = X (T). Therefore X (T) = |x( S ) im P lies 
x{T f ) = |x(Sa,a )- In this case, the identification of the two triangles transforms the 
handlebody H-^ into 7^^ , and the triangles A = Aq is now in the interior of n£ 
By construction, the triangulation T" satisfies H^ Ue aa (-T') ~ M. 

On the other hand, when A is not in the neighborhood of Aq, i.e. E(A,Aq) = 
and V(A, Ao) = 0, 1, the removal of the triangle A = Ao from T creates a hole in T. 
The Euler number of the boundary surface then decreases by two, x(^a,a ) = x(^) ~~ 2. 
On the other hand, the numbers of the edges and the vertices of T' remain the same as 
those of T. Therefore x{T') = x{T) — 1 = |x(^a,a )- The handlebody n^ acquires a 
new handle after the removal of A = Ao, and Ao becomes a meridian disc of the handle. 
Therefore ^( qU ^ a o» u Ea . (-T) ~ M. 

As in the case of eq. ([32]) , we can estimate the contribution of these triangulations to 
Z n3jni (£; a)\Mi taking the symmetry factor C(E)/C(T) into account, and obtain 

C(E Ao ) • 22 ( — - j = 22 Z r t3 ,r tl -E(A,A )( S A,A ;a)|M 

T \ \ / / A.TPfA A _ \_£n 



A:E(A,A )^0 

+ Yl Z n3,ni(^A,Ao',OLU{dA })\ M . 

A:E(A,A )=0 



(33) 



Combining this with eq. (p2|), Z n3jni (E; ol)\m is expressed as 



^n 3 ,ni(S;a) — Z n3 _i jni (E- o ;a) + ^ ^n 3 ,n 1 -E(A,Ao)(SA,A o ; «) 

A:£(A,A )/0 

+ J2 Z n3 , ni {Z AAo ;aU{dA }). 

A:.E(A,A o )=0 



(34) 
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Therefore (T,^ a ')i M solves the Schwinger-Dyson equation, for any choice of M. 



5-2. Triangulations of Handlebodies 

We have seen that the Schwinger-Dyson equation (|34]) holds for any closed orientable 
manifold M. Especially when M = S 3 , there is a nice characterization of triangulations T 
counted in Z n3)ni (E; a)\ S 3- 

If we remove a handlebody from S 3 , the rest is also a handlebody which is described 
as follows. Let us first consider the case when E is isomorphic to a torus. In this case, 
it is convenient to choose a canonical homology basis {a, /3} such that a is the meridian 
of the handlebody (solid torus) H^ and that a and (3 intersects once with each other 
transversely. The cycle (3 is called the longitude of the handlebody. The choice of the 
longitude is unique upto a, i.e. (3 + not is also the longitude of Hj? . Three-dimensional 
manifolds with genus-one Heegaard splittings are classified as 

?4 Q) ue (-?4 nQ+/3) ) ~ s 3 



n^ u E (-n { y a) ) -s^s 2 



?4 a) U E (- H ^ +m0) ) ~ L m , n (ro#0,l) 



(«), 



where L mn is the lense space. Therefore, if H^* U^ (— T) ~ S 3 , T must be a triangulation 
of a handlebody whose meridian (longitude) is the longitude (meridian) of TC 



(a) 



This is also the case for h(E) > 2. Corresponding to the meridians {oti}^ of 
the handlebody Ti^ , we choose simple homology cycles {f3i} on E such that a; and /3; 
(i = 1, ..., /i(E)) make a canonical homology basis, i.e. they have the following intersection 
properties. 

CKi • ttj = , /3i ■ /3j = 

As in the case of torus, the cycles {A} are called the longitudes of Ttj? . It is known in 

:3 



general that the Heegaard splitting of S is unique upto the stable equivalence [p8| , and 
T satisfying 7i^ Ue (—7") should be a triangulation of a handlebody whose meridians 
(longitudes) are the longitudes (meridians) of H^ . 

The sum over triangulations of the handlebody given by (E( a ))|g3 is closely related 
to the sum over triangulations of S 3 . Before explaining the relation, it is instructive to 
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examine the amplitude (£(")) given by eq.(|27|) and its relation to Z q . Let us consider the 
case when £ is homeomorphic to S 2 and is given by a single tetrahedron. For the single 
tetrahedron, the Har tie- Hawking type state $ is given by the 6j-symbol and ri2 = 4. The 
boundary operator A q J2 j^U)F(j; M) for the tetrahedron is then equal to the quartic 
term in the action S(M) upto the factor j^gK q . Therefore, as in the case of the matrix 
model ( |lOl) , the partition function Z q and the amplitude of the tetrahedron are related as 

1 dZ 



{tetrahedron) = 12h. q 



rp \ 



On the other hand, (tetrahedron) is also expanded as a sum over simplicial complexes 
T bounded by the tetrahedron as in eq. (^2|)- Thus the number of closed simplicial 
complexes T weighted with A^ I q (T) and the symmetry factor 1/C(T) is equal to 12n% 
times the number of open simplicial complexes T such that dT = (tetrahedron) weighted 
with A q I q ((tetrahedron) U (— T)) and the symmetry factor (12/C(T)). 

The similar relation holds between (£(")) 153 and the sum over triangulations of S 3 . 
When E is a single tetrahedron, the amplitude is a sum over triangulations of a three- 
dimensional ball B 3 whose boundary is the tetrahedron. 

(tetrahedron) \s 3 = \_\ Z n3 ^ ni {tetrahedron)\s3 ■ g na u ni 

n 3 ,ni 

Here Z n3}ni (tetrahedron) \$ 3 is the number of triangulations of the interior of the tetra- 
hedron with ns tetrahedra and n\ edges weighted with the symmetry factor (12/C(T)). 
The factor 12 is the order of the symmetry of the tetrahedron. We can show that this is 
related to Z n3jUl (S 3 ) as 

Z n3 -i yni (tetrahedron)\ S 3 = 12n 3 Z n3yni (S 3 ). (35) 

Any triangulation of B 3 is obtained by removing a tetrahedron from a triangulation of S 3 , 
and conversely any triangulation of S 3 is given by gluing a tetrahedron to a triangulation of 
B 3 . For each triangulation T of S 3 , there are n%(T) distinct triangulations of B 3 as far as 
T does not have a symmetry. On the other hand, if T has a non-trivial symmetry, different 
choices of tetrahedra in T ~ S 3 related by the symmetry give the same triangulation of 
B 3 . The symmetry factors in Z n3jni (tetrahedron) and Z n3jni (S 3 ) compensate for the 
overcounting due to the symmetries, and eq. ([J5|) holds. 
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6. Factorization 

In the last section, we have shown that (E^)]^ gives a solution to the Schwinger- 
Dyson equation when E(A, Aq) + 1 > V(A, Aq) for any A on £ and when the surface 
£ in the equation is not degenerate. In this section, we will define a set of factorization 
conditions on (H^ a ')\s 3 so that the Schwinger-Dyson equation becomes applicable even 
when it involves degenerate surfaces. We will then show that the Schwinger-Dyson equation 
combined with the factorization conditions has a unique solution. Thus they characterize 
the sum over triangulations T such that Ti^ Us (— T) ~ S 3 . 

Let us first consider a case when the surface £ is degenerate. When the pinching 
cycle is homo logically trivial on £, the surface is actually separated into two components 
Ei and £2 touching to each other at the vertex. In this case, T must also consists of 
two disconnected components T\ and T 2 such that <9Tj = £j (i = 1,2). The Schwinger- 
Dyson equation remains applicable in this case if we supplement a condition that A in 
the summation in the right-hand side of eq. (^4|) is restricted to be those on the same 
boundary component as Ao does. Especially when we can choose the meridians of the 
handlebody Hj? such that 

{«i, ..., «/,(£)} = {a^\ ..., o^)} U {a{\ ..., ajf^} (36) 



as 



where a\ G £1 and a\ G £2, the amplitude (£( a ))is3 factorizes 



(^ ) )| S 3 = (£^ (1) )) |S 3(£^ 2) )) |S 3 



In this case, it is obvious that A G £j should be imposed in eq. (|34| ) when Ao G £i 
(i=l,2). 

It is instructive to compare this behavior of the amplitude with that of (£W) = 
(^$( Q ) 5 \[f E ). When A and Aq wrap around the homologically trivial cycle 7 and when 
the longitudes {a.i\ are separated as in eq. (|36|), «a,a defined by ( |25|) gives a map from 
V^ phys) into V£ h ys) ® V£ h ys) as 

^[$ (a) ]A,A = $ (a<1>) $ (a<2>) . 
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The contribution of the degenerate surface £a,a to the Schwinger-Dyson equation is then 

(4 a) Ao ) 

= (Ei Ql) E^ 2) ) 

_ A |n 2 (Si)+/i(Ei) A |n 2 (S2)+/i(S2) 

x y: * (a(i)) (j (i) )$ (Q(2)) (j (2) )^ 1+S2 (j (i) ,j (2) ) 

0«},0' (2) } 
(« (1) )\/v(« (3) )\ , a V ^(Si)C(S 2 ) ^ n3(T)ffni(T) 



T:9T=Ei+E 2 ^ 



In the first term in the right-hand side, the sum is over simplicial complexes consisting of 
two disconnected components 7\ and T 2 each of which is attached to Ei and £ 2 respectively, 
and the second term is a sum over connected simplicial complexes. We see that the second 
term is multiplied by the extra- factor of A q . This is interpreted as follows. If T remains 
connected after pinching around the cycle 7, T should contain a homotopy cycle linking 
with 7. The extra-factor A q is associated with this homotopy cycle. The factorization 
condition may be regarded as the q — > limit of this equation for (£( a )) where the second 
term is to be suppressed by the factor of A q . 

When the surface E is degenerate around a homologically non-trivial cycle, the sum 
over A in the Schwinger-Dyson equation runs over all triangles on Ea as f ar as ^0 is n °t 
on the degeneration point. When one of the vertices of Ao is on the degeneration point, 
we must supplement the following rule: The sum over triangles in the Schwinger-Dyson 
equation should not include those in the neighborhood of Ao and on the other side of the 
degeneration point across the degenerating cycle. The reason for this is the following. If 
we take a triangle A' which shares the degeneration point with Aq and which is on the 
other side of the point, the identification of A' and Aq does not change the Euler number 
of E. The number of triangles in T, on the other hand, decreases by one if A' and Aq were 
identified in T and removed from T. Thus the equality x(T) = |x(E) cannot be sustained. 
On the other hand, if we take A to be other than those triangles, there is a triangulation 
T such that Ttj? U (— T) ~ S 3 and A and Aq are identified. Thus the Schwinger-Dyson 
equation holds with this supplementary condition. 

So far, we have examined the situation when E is already degenerate. When 
V(A, Ao) > -E^A, Ao) + 2, the triangles A and Ao wrap around a cycle 7 on E and 
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the surface will develop additional degeneration if we identify A and Aq and remove them. 
If it is possible to construct a triangulation T in which the two triangles A and Ao are 
attached to each other, T should be counted in (E( a ))j S 3. The Schwinger-Dyson equation 
( p6|) derived for non-degenerate surfaces becomes applicable in this case if we take into 
account the contribution of the degenerate surface Sa,a - On the other hand, if there is 
no such triangulation in which A is identified with A , the surface Ea,a should not be 
included the Schwinger-Dyson equation. 

When the cycle 7 is homologically trivial, it is always possible to construct a triangu- 
lation T in which A and Ao are identified. Thus the term (£^A ) should be included in 
the Schwinger-Dyson equation. In this case, the surface Ea,a is pinched around at the 
common vertex of A and Ao, and the Schwinger-Dyson equation can be iterated with the 
supplementary conditions introduced in the above. 

On the other hand, when the cycle 7 is homologically non-trivial on E, the triangles 
A and Ao wrap around a handle on E. In this case, the handle is pinched on the surface 
Sa,a - If there is a triangulation T in which the triangles A and Ao are identified with 
each other, there must be a disc in T whose boundary is the cycle 7. After removing the 
triangle A = Ao from T, we obtain T' which gives a triangulation of a handlebody whose 
boundary is Ea,a - The surface Ea,a is pinched at the common vertex of A and Ao. 
The complement of T' in S 3 should also be a handlebody 7i^ which is pinched around 
the cycle 7. The process TC^ U (— T) — ► H^ U (— T') then gives the reduction of the 
Heegaard splitting of S 3 which is the inverse of the stabilization ||38|| . For this to be possible, 



there must be a disc in the handlebody 1~C£ whose boundary is on E and intersects once 
with 7 transversely. If this is the case, we can choose the meridians of T~tj? so that a±-j — 1 
and a, -7 = (i = 2, ..., h(T,)). The interior of 7^ E 1 '"' MS) is then homeomorphic to that 
of 7iyv 2 ' MS) , where E' A A is obtained from Ea A n by splitting the common vertex of A 
and Ao into two so that the pinched handle is removed; h(E' A A ) = /i(Ea,a') — 1- In this 
case, a term (E A ° A )m must be included in the Schwinger-Dyson equation. Conversely, if 
it is not possible to find such a disc in 7i E , there is no triangulation in which A and Aq are 
identified with each other. In this case, there should be no term in the Schwinger-Dyson 
equation corresponding to the surface Ea,a - 

With these supplementary conditions, we can apply the Schwinger-Dyson equation 
even when E is degenerate or when there is a triangle A such that V(A, Ao) > E(A, Ao)+2. 
As in the case of the matrix model, we can show that the Schwinger-Dyson equation com- 
bined with these conditions can be solved inductively. The Schwinger-Dyson equation 
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expressed as in eq.([J3|) is a recurrence relation for Z n3jni (a; E) along the bibliographical 
ordering: (77,3, E) y (n 3 , E') if 77-3 > n 3 or 77,3 = n 3 and 77.2(E) > tt^E'). Because of the fac- 
torization conditions, we can trace back the Schwinger-Dyson equation to the expectation 
value of the null surface (1) which is normalized to be one. Thus the Schwinger-Dyson 
equation combined with the factorization conditions has a unique solution. The solution 
gives the sum over triangulations of handlebodies, and it is related to Z n3 , ni (S 3 ) as in eq. 

(PD- 
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